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Abstract 

In this paper we review, and extend to the non-isothermal case, some results con¬ 
cerning the application of the maximum entropy closure technique to the derivation of 
hydrodynamic equations for particles with spin-orbit interaction and Fermi-Dirac statis¬ 
tics. In the second part of the paper we treat in more details the case of electrons on a 
graphene sheet and investigate various asymptotic regimes. 

1 Introduction 

This paper is devoted to present some results on the derivation of hydrodynamic equa¬ 
tions describing electrons subject to spin-orbit-like interactions. Systems of this kind, of 
particular interest for applications to microelectronics, include electrons undergoing the 
so-called Rashba effect [T], the Kane’s two-band K-P model [5] and electrons in single¬ 
layer graphene [3]. The diffusive and hydrodynamic descriptions of such systems are 
extensively treated in Refs. [1 El El Elia Ellin]- Here, we summarize the results contained 
in Refs. Enaiann], concerning the hydrodynamic description, and extend them to the 
non-isothermal case. 

In comparison with kinetic models, the advantages of fluid models for applications 
are evident. In fact, from the numerical point of view, a system of PDFs for a set of 
macroscopic quantities is much more desirable than a single equation for a density in 
phase-space, where also the components of momentum are independent variables. More¬ 
over, from the point of view of mathematical modeling, they offer more flexibility, as 
various kind of boundary conditions and coupling terms (e.g. with a self-consistent po¬ 
tential or with various types of scattering mechanisms) can be very naturally embodied 
in the mathematical model. 

On the other hand, the derivation of fluid equations for the systems under consider¬ 
ation, which possess spinorial degrees of freedom and non-parabolic dispersion relations 
(energy bands), is far from being a trivial extension of the techniques employed for stan¬ 
dard (i.e. scalar and parabolic) particles. The best strategy to obtain hydrodynamic (or, 
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more in general, fluid-dynamic) equations in this case, seems to be their systematic deriva¬ 
tion from an underlying kinetic description by means of the Maximum Entropy Principle 
(MEP) and its quantum extensions [11] |T2 [131 [14] . The MEP basically stipulates that the 
microscopic (kinetic) state of the system is the most probable among all states sharing the 
same macroscopic moments of interest, providing therefore a formal closure of the system 
of moment equations. This is a very general principle which finds a variety of applications 
to different fields, ranging from statistical mechanics to signal theory [15] . For quantum 
systems it can be used in combination with the quantum kinetic framework provided by 
the phase-space formulation of quantum mechanics due to Wigner [lain]. 

In the present work the Wigner formalism is used “semiclassically”, which means that 
some quantum features are retained (namely, the peculiar energy-band dispersion relations 
and the Fermi-Dirac statistics) while others are neglected (namely, the quantum coherence 
between different bands). Consequently, the obtained hydrodynamic description misses 
some interesting physics when quantum interference between bands becomes important 
(e.g. close to abrupt potential variations [Tang). Nevertheless, the derived equations 
possess an interesting mathematical structure and reveal some interesting physics, still 
occurring in absence of such interference phenomena (see, in particular. Section [S]). 

The paper is organized as follows. In Sectionthe kinetic-level formalism, based on a 
semiclassical Wigner description, is introduced for a fairly general spin-orbit Hamiltonian 
that includes all cases of interest. In Sec. ia we write the moment equations for density, 
velocity and energy, and perform their formal closure by means of the MEP. Then, the 
second part of the paper is focused on the case of graphene. In Sec. H] the general theory 
exposed in the first part is specialized for the Dirac-like Hamiltonian describing electrons 
on a single-layer graphene sheet. In Sec. |S]we obtain the asymptotic form of the hydro- 
dynamic equations derived in Sec. |4] in some physically relevant limits (namely, the high 
temperature, zero temperature, collimation and diffusive limits). Finally, Sec.|n]is devoted 
to conclusions and perspectives. 

2 Phase-space description of spin-orbit particles 

Let us consider a spin-orbit Hamiltonian of the form 

^(x,P) = [/lo(p) + l^(x)]cro + h(p) • (7, (1) 

where x S K^, p S K^, cr = (cti, tT 2 , 0 - 3 ), h = (/ii, / 12 , hs) and 



Moreover, the dot product is defined as h • cr = hiai + h 2 <J 2 + 

Hamiltonians of this kind describe various systems of great interest in solid-state 
physics. The first example is a 2-dimensional electron gas confined in an asymmetric 
potential well, which is subject to the Rashba spin-orbit interaction mm- In this case: 

1 2 

d = 2, ho{p) =-;—\p\ , h(p) = apxe;,, 

2m* 
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where p = {px,Py), is the electron effective mass, is the normal direction to the 
well and a is the Rashba constant. 

Another example is the two-band K-P model [aiail0l[20]; in this case: 

d = 3, /lo(p) = ^|p|^ h(p) = ^0, ^K-p, , 

where m is the electron (bare) mass, Eg is the band-gap and K is the matrix element of 
the gradient operator between conduction and valence Bloch functions. 

The last example is that of electrons on a single-layer graphene sheet [a [3 [n], in 
which case: 

d = 2, ho{p) = 0, h(p) = cp. 


This case will be considered in more details in the second part of the paper. 

We remark that the variable p has to be interpreted as the crystal pseudo-momentum, 
rather than the ordinary momentum. The interpretation of the vector variable a depends 
on the cases: it is (proportional to) the spin vector in the case of Rashba Hamiltonian, 
while it is a pseudo-spin in the other two examples 013]. For graphene, in particular, the 
pseudo-spin is related to the decomposition of the honeycomb lattice into two inequivalent 
sublattices, which reflects the presence of two carbon atoms in the fundamental cell of the 
lattice 0121]. 

The main semiclassical quantities associated with o are: 


1 . the two energy bands 

^^±(p) = ^o(p) ± |h(p)| 

(i.e., the eigenvalues of H with ]/ = 0); 

2. the projectors on the eigenspaces corresponding to E±{p), 


P±{P) = 2 (^ 0 ±J^(P) -cr), 


where 


v{p) = 


h(p) 

|h(p)| 


is the pseudo-momentum direction; 

3. the semiclassical velocities 

v±(p) = VpF;±(p); 

4. the effective-mass tensor m 


( 2 ) 

(3) 

(4) 

(5) 


M^^(p) = Vp (g) v±(p) = Vp (g) Vp£;±(p). 


( 6 ) 


Note, in particular, that the eigenvalues of the projector P± are 1 and 0, corresponding 
to whether or not the electron energy belongs to the upper/lower energy band. Hence, 
the expected value of P± can be interpreted as the fraction of electrons belonging to the 
upper/lower band (see below). 

The phase-space description of a statistical population of electrons with Hamiltonian 
O is provided by the Wigner matrix [illllllTlIlS] 

3 

F(x,p,t) = y]/fc(x,p,f)CTfc, (7) 

/c=0 


3 



which is the Wigner transform, 

/fc(x,p,t)= [ pfe (x+^,x- 

jRd \ Z Z/ 

of the spinoriai density matrix 

3 

p(x, y, 0 = X! 

fc =0 

Such representation of a quantum mixed-state has the fundamentai property that the 
expected vaiue of an observabie with symboi A = X]fc=o p)o’fe is given by the ciassicai- 
iooking formuia 

3 

Ef[A]= f Tr(FA)dxdp = ^ f V afe(x, p)/^(x, p, t) dxdp. (8) 

dR2<i (27rfi,)“ Jg2<i 

By appiying Eq. ([5]) to the band projectors ^’±(p) we obtain 

]Ef[-P±] = . [ (/o ± 1 / • f) dxdp 

{2TThf J^2d 

and it is therefore naturai to interpret the functions 

/± = /o ± ■ f (9) 

as the phase-space densities of eiectrons having energies, respectiveiy, in the upper and 
iower band. 

Let us now consider the foiiowing hydrodynamic moments of eiectrons in the two 
bands: 

n± = (/±), (density), 

n±u± = (v±/±) (veiocity), (10) 

n±e± = {E± f±) (energy), 

(see aiso Ref. [M] where additionai moments are considered). Here we have introduced 
the shorthand 

if) (x, t) = ^ f /(x, p, t) dp. 

i^TThf 

The (semiciassicai) dynamics of the Wigner matrix (jT]) is provided by the Wigner equations 
for the Hamiitonian (HD m, 

3 

{dt + Vp/iQ • Vx — VxH • Vp) fo -f Vp/ifc ■ Vx/fe = 0, 

< ( 11 ) 

(dt + Vp/io • Vx - VxH • Vp) f, + Vph, ■ Vx/o = -^ (h x f)i, 

with i = 1,2, 3. From (ITT]) , the foliowing equations for the band-Wigner functions /+ and 
/_ (see definition ®) are readily obtained: 

{dt + v± • Vx - VxE • Vp) /± = -Vx • f_i ± F • (VxH • Vp)fx, (12) 

where the terms containing 

fj_ := (z/ X f) X F 

are responsible for quantum interference between the two bands [ai23]- 
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3 Maximum entropy closure 


In order to obtain from (El a closed system of equations for the moments (nni , we 
assume that the system is in a state of maximum entropy, according to the so-called 
Maximum Entropy Principle (MEP) [TTJ [T^l [131 [IH HI] which in the present case reads 
as follows: 


MEP F™® is the most probable microscopic state with the observed macroscopic moments 
n±, u± and e±. 

Hence, we search for a Wigner matrix F™® that maximizes the total entropy 


£{F) = - / Tr{s(F)}(x,p)dpdx 

(ZTTfi.) J-g2d 

among all matrices F = J2t=o fk<^k, such that 0 < F < 1 and 





V± 

W) 

f±) = n± 

u± 

WJ 


(13) 


(14) 


In (1131) . ks is the Boltzmann constant, Tr is the matrix trace and 


s(a:) = X log X-I-(1 — x) log(l — x), 0 < x < 1, (15) 

is (minus) the Fermi-Dirac entropy function. The condition 0 < F < 1 ensures that s(F) 
is a well-defined matrix. 

It can be proven [5] that 


/r = (s')-' (v± • B± + - C±F±) 

exp(C±F± - v± • B± - H±) -I- 1’ 

where A±, B± = {Bi,... ,Bd)± and C± are Lagrange multipliers (functions of x and t), 
and, moreover, 

= 0. (17) 

Thus, the (semiclassical) MEP state corresponds to two local Fermi-Dirac distributions in 
the two energy bands. In particular, Eq. (113 implies that, in such state, the interference 
terms vanish and, therefore, the two bands are decoupled (unless additional coupling 
mechanisms are considered El Hi HI). Hence, from now on, we shall treat the two bands 
separately and, in order to simplify notations, the ± labels will be suppressed (except 
where a distinction between quantities taking different forms in the two bands, such as 
E± or v±, is necessary). 

Using m and (HB in (fT^ (and suppressing the ± labels, as it was just explained) 
yields 

(dt + v± • Vx - VxU • Vp) = 0. (18) 
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By taking the moments (•), (v±-) and {E±-) of both sides of Eq. (fT5|) . and recalling the 
definitions and we obtain the moment equations 

{ dtn + dj(nuj) = 0, 

dt{nui) + djP^ + QfjdjV = 0 , ( 19 ) 

dt{ne) + djS^ + nujdjV = 0 

where di = djdxi and 

= {vtvffn, 

Qti = m ri, ( 20 ) 

sf = {E±vf r^). 

Thanks to the MEP, the moment system (HU) is implicitly closed by the constraints 
dMD, linking the Lagrange multipliers (A, B, C) to the moments (n, u, e) thus allowing (in 
principle) to think to as being parametrized by (n, u, e) and, consequently, the extra 
moments (P^, Q^, S^) as functions of the unknowns (n, u, e). 

Following Levermore we can express the moments of the MEP state as the 
derivatives with respect to the Lagrange multipliers of the “density potential” e*, which 
is defined as the Legendre transform of the entropy density 


e={s{ri) 


where s is given by (fTSl) and by (ITBl) . It is not difficult to show that 


£* = -(log(l-/-*=)) 


and that the constraint equations HI may be rewritten as 


de* 


dA 


de* 

OB, 


= nui, 


de* 
' dC 


= ne, 


( 21 ) 


where i = Levermore’s theory, moreover, ensures that system HI, with the 

closure relations (uni) is hyperbolic and, therefore, it is at least locally well-posed (see also 
Ref. H). 


4 The case of graphene 

We now specialize the formalism introduced so far to the case of a population of electrons 
on a single-layer graphene sheet. Such electrons, in the proximity of a Dirac point in 
pseudo-momentum space [3], are described by the Hamiltonian o with 

d = 2, ho{p) = 0, h(p) = cp, 

(where c ~ 10®m/s is the Fermi velocity), which corresponds to a Dirac-like Hamiltonian 
for relativistic, massless particles. We remark that this is an approximation which is valid 
only in the proximity of a Dirac point for an infinite, ideal and un-doped system (see Ref. 
m and references therein). In this case the energy bands are the Dirac cones 

E±(p)=±c|p|, (22) 
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and the eigenprojections are given by 

P±ijp) = \{.^o±v{jp)-(j), (23) 

where 

Kp) = (24) 

Moreover, the semiclassical velocities are 

v±(p) = ±-^ ± c: 2 (p), (25) 

implying that electrons travel with the constant speed c and direction v, and the effective- 
mass tensor is 

(p) = (p) ® (p) (26) 

where 

= {-V2,Vi). 

Since the lower band is unbounded from below, we have to change a little the theory 
developed in the previous sections and describe the lower-band population in terms of 
electron vacancies, i.e. holes. This is achieved by means of the substitution 


/_(x,p,t) I— >1- /_(x,-p,t), 

which brings the transport equation, Eq. m , into 

(ai+ci/-V xTVxE-V p)/”*^ =0. (27) 

Note that the only difference between electrons and holes is the charge sign. Moreover, 
the MEP-states for electrons and holes have now the form 

^ exp(C|p|-Kp)-B-A) + l’ 

in fact, both upper-cone electrons and lower-cone holes have positive energies 

E{p) = c|p| (29) 

(note that in (E7)) the Fermi velocity c has been absorbed in the Lagrange multiplier C). 
Moreover, we slightly change the definition of u to be the average direction 

nu = (j 2 /), 0 < |u| < 1, 


which differs from average velocity just for the constant factor c. The inequality |u| < 1 
is an obvious consequence of the fact that u is an average of directions. 

The moment equations (USD, in the specific case of graphene, read as follows: 

{ dtn + cdj{nuj) = 0, 

dt{nui) + c^jP^j ± Q^j^jV = 0, (30) 

dt{ne) + cdjSj ± cnujdjV = 0, 
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where the higher-order moments Pij, Qij and Sj take the form 




1 


Sj = {cpjfn- 


( 31 ) 


We now intend to find an (as much as possible) explicit expression for the dependence 
of the Lagrange multipliers yl, B = (Bi, B 2 ) and C in terms of the moments n, u = (rti, U 2 ) 
and e, as resulting from the constraint equations 

{ri = n. {c\v\n) = ne. (32) 

By expressing the integrals over p S in polar coordinates, we obtain the expressions 

Io(A|B|) 


in = 


{I'ri = 


{c\p\fn = 


27rfi2C'2 ’ 

pliAm B 

27r/i2C'2 |B|’ 

ciiiAm 


(33) 




where 


1 r 

I%{x,y) = — / cos{N0)(j)s{x+ ycos0)d0, 
^ Jo 

and (ps is the Fermi integral of order s > 0: 

1 /■“ P- ' 


(34) 


It is now convenient to put 
B = IBL 


AA) = 


c = 


r(s) 


+ 1 


• dt. 


riT = 


fc|r 2 


1 


keT' 27r/i2c2 

so that the previous expressions can be rewritten as 

{ri=nTXl{AB) 

An = '^PliAB)B, 
{c\p\n=2nTkBTll{A,B). 


27rfi2C'2 ’ 


(35) 


(36) 


We remark that the new Lagrange multiplier T has the physical meaning of the electron 
gas temperature. From (1551) and the constraint equations (1551) . we obtain that B has 
the same direction as u and that (n, |u|, e) are related to the scalar Lagrange multipliers 
{AB,T) by 

ll{A,B)nT = n, 


II{A,B) 

PliAB) 

PliAB) 

Tl{A,B) 


= u 


2 fcijT = e. 


( 37 ) 












Similarly to what is found in Ref. [7], we obtain the following expressions of the higher- 
order moments m in terms of n, u, T and the functions = I%{A, B): 




(38) 



2kBTn l\ 


where = (—U 2 ,iti). 

5 Asymptotic regimes 

The expressions (1381) of Qij and Sj are still not explicit, as functions of n and u. In 
fact, these expression depend, through the functions I%{A, B), on the two scalar Lagrange 
multipliers A and B, which are related to n and |u| via the relations (1371) . In Ref. [7] it 
has been proven that the correspondence between {A,B) and (n, |u|) is 1-1 but, as far 
as we know, it is not possible to give an explicit, analytic, expression of the former as 
functions of the latter. 

However, we can say more in some particular regimes of physical interest. Such regimes 
correspond to different asymptotic regions [7] in the half plane (H, B) S Rx [0, oo), namely: 

1. the asymptotic region + B"^ ^ oo with A < —B (i.e. {A, B) below the “critical 
line” A + B = Q), corresponds to a regime of high temperatures, where the Fermi- 
Dirac distribution is well approximated by a Maxwell-Boltzmann distribution; 

2. the asymptotic region A^ + B^ ^ oo with A > —B corresponds to the limit T —0, 
in which case we speak of “degenerate fermion gas”; 

3. the asymptotic region -|- —>■ oo with A ^ B (i.e. {A, B) approaches the critical 

line A + B — 0) corresponds to a “collimation regime”, |u| —>■ 1, where the velocities 
of the electrons are all aligned along a ((x, ^(-dependent) direction in the p-space 
(the direction determined by u); there are two types of collimation, depending on 
whether the critical line is approached from below (Maxwell-Boltzmann collimation) 
or from above (degenerate gas collimation); 

4. opposite to the collimation limit, the asymptotic region P —^ 0 corresponds to the 
diffusive limit |u| —> 0, where the velocities are randomly spread over all directions. 

The asymptotic analysis of equations (15(11) in these regimes is based on the following 
result, which has been proven in Ref. [7]. 

Theorem. The functions have the following asymptotic behavior: 

1. in the Maxwell-Boltzmann limit, -f P^ —>■ oo, with A < —B, 


11,{A,B) ^ In{B), 


( 39 ) 


where In are the modihed Bessel functions of the first kind; 
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2. in the degenerate gas limit, ^ oo, with A > —B, 


where 




7rr(s + 1 ) 


fCiA,B) 

cos{N0){A ■ 


B cos 0yd0, 


C(AS) = 5? [cos-1 (-|)] 


arccos (—^) , if —B < A < B, 
TT, if > B. 


(40) 


(41) 


5.1 Maxwell-Boltzmann regime 

The Maxwell-Boltzmann regime is the limit for large T and corresponds to + B'^ —)> oo 
with A < —B in the {A,B) half plane. In this case, we can use the approximation (l39l) . 
Note, in particular, that in such limit the functions become factorized and independent 
on the index s. Then, the constraint equations (I37p become 


e^Io{B)nT = n. 


h{B) 

Io{B) 


2fcijT = e. 


(42) 


and it can be shown that the the MEP-state (1^51) is well approximated by the Maxwellian- 
like distribution 


jme _ 


riT Io{B) 


exp 


.^|p|+B^(p). " 

KbT u 


where 



Moreover, we get the explicit form of Pij, Qij and Sj: 

[^(|u|)MiMj + (l-^(|u|))u,^u7^] , 

|u| 

'3d = [^(|u|)u,^uj--b(l-X(|u|))M,Uj] , 

e |u| 


(43) 

(44) 


(45) 


Sj = neuj 


where 


^(|u|) 


Io{B)+l2iB) 

2/o(B) 


and B is given by (HT)) . 

By playing a little with the asymptotic expansions of the modified Bessel functions /„ 
we obtain the asymptotic behavior of X(|u|) in the diffusive limit: 


^(|u|) = i + i|u|Ve)(|u|'‘), as|u|^0 , 

and in the collimation limit: 


(46) 


X(|u|) = 1-2(1 -|u|)2 + 0((l-|u|)3), as|u|^l. (47) 
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Substituting Sj = neuj in the third of the moment equations (1501) yields, after a little 
algebra, 

dte + cujdjc ± cUjdjV = 0. (48) 

Thus, the isothermal case (e = 2 A:bT constant) is only compatible with UjdjV = 0, i.e. 
the component of the force field parallel to the velocity field must vanish. In this case, 
the pseudo-momentum balance equation reduces to 


dt{nui) + cdj 


nX(|u|) UiUj 


U 


n{l- X{\u\))u:tuf \ 

|u|^ ) 


± = 0. (49) 


5.2 Degenerate gas regime 

The degenerate gas regime is the limit for T —>■ 0 and corresponds to ^ oo with 

A > —B, in the (A, B) half plane. In this case, we can use the approximation (1401) "(l4T |) . 
It is convenient to put 

A = Rcos4’, B = Rsinip, 

and rewrite (H 0 |) as follows: 


Stt 


where 


and 


= 


7rr(s -I- 1) Jq 


rcm 

cos{N0){cos'4) + cos OydO 


C(')/')= 5R [cos ^(—coti/')] = 

The asymptotic form of the constraint equations (|57)) is now 


TT 

4 


if 0 < -0 < —. 

4 


' 5 

(50) 


(51) 

Stt 


T’ 

TT 

(52) 


J'q (■ 0 )i?^n 7 ’ = n, 




= U 


ksTR = e. 


(53) 


Note that: 

1 . |u| only depends on ij) and we can write 

(|u|); (54) 


2. recalling (1551) . from the first of the above equations we have R ^ 1/T and, then, the 
third equation shows that e remains positive even though T —^ 0. 


II 








From the above considerations it is readily seen that, in the limit T —>■ 0, the MEP-state 
(I^Rll takes the typical degenerate Fermi-Dirac form 


jme _ Q 


12'Kh?n 

HW) 


^^(p) 


U 


■ sin Ip + cos tp 


(55) 


where 9 denotes the Heaviside function and '0(u) is given by (IMll . Using dSS]), (1^ and 
we obtain the following expressions of Pij and Qij and Sj for a degenerate electron 

gas: 


P — 

JTij — 

Qij — 


[y(|u|)uiUj + (i-y(|u|))M,^M^^], 


[Z{\u\)u^Uj + Z±{\u\)uiuf] , 


hy/Tr\u 


Sj = lU(|u|) 


neuj 


(56) 


where 


YM)-- 

^±(|u|) 


2F^{iP) ’ 


Z{H) 


nii’) - w 

2./2T^ 


lU(|u|) 




and Ip = ■i/'(|u|) is given by Eq. (15^ . 

By using the techniques developed in Ref. [7] , it is not difficult to calculate the asymp¬ 
totic behavior of the functions y(|u|), .Z(|u|), Zj_(|u|) and lU(|u|) in the two limits |u| —?> 0 
(diffusion) and |u| —>■ 1 (collimation). 

For |u| —>■ 0 we obtain: 


y(|u|) = i + i|u|^ + o(|uh, 
Z(|u|) = i-i|u|VO(|uh, 
^u(|u|) = i-i|u|VO(|uh, 

lU(|u|) = ^|u|+0(|u|3). 


(57) 


For |u| —>■ 1 we obtain 


y(|u|) = l-2 (l-|u|) + 0((l-|u|)2), 

Z(|u|) = -^(l-|u|)t-fO((l-|u|)?), 

Z^(|u|) = (1 - |u|)^ + 0((1 - |u|)t), 

lU(|u|) = l-I(l-|u|) + 0((l-|u|)2). 
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5.3 Collimation regime 


The collimation limit corresponds to the absence of spread in the particle directions, 
i.e. to |u| —y 1. It can be shown [315] that this limit is equivalent to ^ oo 

with A/B —>■ —1. However, there is a completely different behavior when the critical 
line A = —B is approached from below (Maxwell-Boltzmann collimation) of from above 
(degenerate gas collimation). 

The first case corresponds to taking the limit H —>■ oo in the “Maxwellian” distribution 
(15^ . which produces a delta in the angle between p and u, namely 


jme 


27rn 

-exp 

riT 
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(59) 


Moreover, since -^(|u|) —>• 1 as |u| —>• 1 (see Eq. (03)), from Eq. (1^ we obtain 

2 c I I 

Pij —y TlUiUj,, Qij —y TlUj^ Uj 

and the pseudo-momentum balance equation reduces to 

2 c 

dt{nui) -\- cdj{nUiUj) ± — nu^uj~djV = 0 . 

By using the continuity equation dtn + cdj{nuj) the latter can be rewritten as 

2 c 

dtUi + cujdjUi ± —ufufdjV = 0 (60) 

which is decoupled from the continuity equation for n. As pointed out in Refs. [3IH], this 
equation reveals that collimated electrons in graphene have the properties of a geometrical- 
optics system, with “refractive index” 


Af(x) = e 


^ pT|V(x) ^ gT, 


,V(x) 


By also considering the energy balance equation (l48ll . we finally obtain the system 


{ dtUi + cujdjUi ± —djV = 0, 

dte + cujdjC ± cujdjV = 0 . 


(61) 


In order to derive the collimation equations for a degenerate gas, we start from the 
expression (1551) of Pij , Qij and Sj , and use the asymptotic relations (1551) to obtain that 


Pij —>■ nuiUj, Qij —0, Sj —>■ neuj, 

as |u| —1. But then, the hydrodynamic system (1551) degenerates into the decoupled 
system 

dtn + cdj{nuj) = 0, 

dtUi -\- cUjdjUi = 0 , 

dte + cujdje ± cUjdjV = 0 . 

Such a “trivial” asymptotic behavior of collimated degenerate electrons has already been 
pointed out in Ref. |3 in the isothermal case, and is due to the vanishing effective-mass 
tensor Q. 
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5.4 Diffusion regime 


The diffusion regime corresponds to the limit |u| —> 0 of vanishing mean velocity. In order 
to observe the diffusive behavior we have to introduce in (I30p a current-relaxation term 
—nu/r, and to rescale time and velocity as 

t* = Tt, u* = —u. 

T 

In this way we obtain the system 


dt*n cdj(nu*) = 0, 

df inu*) + cdjPij ± QijdjV = -nu* (62) 

rdt* (ne) -I- cdjSj ± CTnu*djV = 0, 


where the terms Qij and Sj, depending only on u/|u| = u*/|u*|, remain unchanged 
except that the Lagrange multipliers must satisfy 




r u 


(63) 


In the diffusive limit r —>■ 0 we obtain the condition X^{A^B) = 0, which is satisfied if 
and only if P = 0 [7]. Since 


I%{A, 0 ) 


MA), ■an = o, 

0, if > I, 


from the first of (1571) with P = 0 we obtain 


(64) 


A = (65) 

and, moreover, 

Tl 

Pij{A, ~ j 

^^(AO) = 0 . 

Letting r —0 in Eq. (1621) yields, therefore, the diffusive system 


dt*n + cdj(nu*) = 0, 
nu* = - {cdjPij ± QijdjV), 


(67) 


with P = P(4,0) and Q = Q{A,0) given by (| 66 )) . that is, in terms of the original time 
variable 


dtn=^d, 




( 68 ) 


It is not difficult to check that the diffusion equation (l68l) take the specific form 



(djU ± , ^ 

V ^ ksT ^ J 


(69) 
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in the Maxwell-Boltzmann limit and 


dtn = ^ dj (^cdj-n ± . ( 70 ) 

in the degenerate gas limit. 

We remark that, owing to the conical dispersion relation (1^ . the drift-diffusion equa¬ 
tions (EHH, (EH) and cni) have a “specular” structure with respect to the drift-diffusion 
equations for Fermions with the usual parabolic dispersion relation [niiMiiis]. Indeed, 
the diffusion coefficient (which is proportional to the variance of the velocity distribution), 
is here independent of the temperature T, because the particles move with constant speed 
c, while it is proportional to T in the parabolic case. On the other hand, the mobility 
coefficient (which is related to the distribution of the second derivative of the energy, i.e. 
to the effective-mass tensor) is here temperature-dependent while in the parabolic case is 
constant. Also the nonlinearity, which in the parabolic case affects the diffusive term, in 
Eqs. (IH71) and cni) is found in the drift term. 

6 Conclusions 

We have presented the systematic derivation from the Maximum Entropy Principle of 
hydrodynamic equations describing a population of electrons subject to spin-orbit inter¬ 
actions. In the second part of the paper we have treated more extensively the case of 
electrons on a single-layer graphene sheet. 

The hydrodynamic equations have the form of a Euler-like system of conservation laws 
for density, n, momentum, u, and energy, e, in each of the two bands (the band indices are 
here omitted). Such system is of hyperbolic character, which ensures its (at least) local 
well-posedness. It is worth to remark that the full nonlinear structure of the MEP-state 
is retained, so that no assumptions of linear response or quasi-isotropic distribution are 
needed. 

The system, in general, is not explicitly closed, i.e. no explicit constitutive relations, 
expressing the higher-order moments Pij, Qij and Sj as functions of n, u and e, can be 
given. However, in the case of graphene and for particular asymptotic regimes (namely, 
the limits of high and zero temperature, the limit of collimated direction and the diffusive 
limit), the closure is fully explicit. 

As already mentioned in the Introduction, our results are not able to capture the 
physics of the system when the semiclassical approximation is not valid, that is when 
the quantum coherence becomes important. This typically happens in presence of rapid 
potential variations, such as potential steps or barriers. In these cases one expects the 
equations derived here to be a good approximation in a “semiclassical region”, far enough 
from the potential steps (constituting instead the “quantum region”). The semiclassical 
regions could be coupled to the quantum ones by means of quantum-classical interface 
conditions, analogous to those developed for standard, i.e. scalar and parabolic, particles 
(see Ref. m and references therein). 


15 



Acknowledgements 

This work has been partially supported by INdAM-GNFM, Progetto Giovani Ricercatori 

2013 Quantum fluid dynamics of identical particles: analytical and numerical study. 

References 

[1] Zutic, L, Fabian, J., Das Sarma, S.: Spintronics: fundamentals and applications. Rev 
Mod Phys 76, 323-410 (2002) 

[2] Kane, E.O.: The k-p method. In: Willardson, R.K., Beer, A.C. (eds.) Physics of 
III-V Gompounds, Semiconductors and Semimetals vol. 1, chap. 3. Academic Press, 
New York (1966) 

[3] Gastro Neto, A.H., Guinea, F., Peres, N.M.R., Novoselov, K.S., Geim, A.K.: The 
electronic properties of graphene. Rev Mod Phys 81, 109-162 (2009) 

[4] Barletti, L., Mehats, F.: Quantum drift-diffusion modeling of spin transport in nanos¬ 
tructures. J Math Phys 51, 053304 (2010) 

[5] Barletti, L., Frosali, G.: Diffusive limit of the two-band k-p model for semiconductors. 
J Stat Phys 139, 280-306 (2010) 

[6] Possanner, S., Negulescu, C.: Diffusion limit of a generalized matrix Boltzmann 
equation for spin-polarized transport. Kinet Relat Mod 4, 1159-1191 (2011) 

[7] Barletti, L.: Hydrodynamic equations for electrons in graphene obtained from the 
maximum entropy principle. J Math Phys 55, 083303 (2014) 

[8] Morandi, O., Barletti, L.: Particle dynamics in graphene: collimated beam limit. 
Journal of Gomputational and Theoretical Transport 43, 1-15 (2014) 

[9] Barletti, L., Frosali, G., Morandi, O.: Kinetic and hydrodynamic models for multi¬ 
band quantum transport in crystals. In: Ehrhardt, M., Koprucki, T. (eds.) Multi- 
Band Effective Mass Approximations: Advanced Mathematical Models and Numer¬ 
ical Techniques. Springer Verlag, Berlin (2014) 

[10] Barletti, L., Borgioli, G., Frosali, G.: Semiclassical hydrodynamics of a quantum 
Kane model for semiconductors. Trudy Instituta Matematiki 11, 11-29 (2014) 

[11] Levermore, C.D.: Moment closure hierarchies for kinetic theories. J Stat Phys 83, 
1021-1065 (1996) 

[12] Degond, P., Ringhofer, G.: Quantum moment hydrodynamics and the entropy prin¬ 
ciple. J Stat Phys 112, 587-628 (2003) 

[13] Trovato, M., Reggiani, L.: Quantum maximum entropy principle for a system of 
identical particles. Phys Rev E 81, 021119 (2010) 

[14] Gamiola, V.D., Romano, V.: Hydrodynamical model for charge transport in 
graphene. J Stat Phys 157, 1114-1137 (2014) 

[15] Wu, N.: The Maximum Entropy Method. Springer Verlag, Berlin (1997) 


16 



[16] Zachos, C.K., Fairlie, D.B., Curtright, T.L. (eds.) Quantum Mechanics in Phase 
Space. An Overview with Selected Papers. World Scientific Publishing, Hackensack 
NJ (2005) 

[17] Barletti, L.: A mathematical introduction to the Wigner formulation of quantum 
mechanics. Boll Unione Mat Ital B 6-B(8), 693-716 (2003) 

[18] Katsnelson, M.L, Novoselov, K.S., Geim, A.K.: Chiral tunnelling and the Klein 
paradox in graphene. Nature Physics 2(9), 620-625 (2006) 

[19] Cheianov, V.V., Fal’ko, V., Altshuler, B.L.: The focusing of electron flow and a 
Veselago lens in graphene. Science 315, 1252-1255 (2007) 

[20] Barletti, L., Ben Abdallah, N.: Quantum transport in crystals: effective-mass theo¬ 
rem and K-P Hamiltonians. Comm Math Phys 307, 567-607 (2011) 

[21] Deretzis, L, La Magna, A.: Origin and impact of sublattice symmetry breaking in 
nitrogen-doped graphene. Phys Rev B 89, 115408 (2014) 

[22] Slonczewski, J.C., Weiss, P.R.: Band structure of graphite. Phys Rev 109, 272-279 
(1958) 

[23] Morandi, O.: Wigner-function formalism applied to the Zener band transition in a 
semiconductor. Phys Rev B 80, 02430 (2009) 

[24] Jiingel, A., Krause, S., Pietra, P.: Diffusive semiconductor moment equations using 
Fermi-Dirac statistics. Z Angew Math Phys 62, 623-639 (2011) 

[25] Barletti, L., Cintolesi, C.: Derivation of isothermal quantum fluid equations with 
Fermi-Dirac and Bose-Einstein statistics. J Stat Phys 148, 353-386 (2012) 

[26] Degond, P., El Ayyadi, A.: A coupled Schrddinger drift-diffusion model for quantum 
semiconductor device simulations. J Comput Phys 181, 222-259 (2002) 


17 



